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Small-Scale Kinematic Dynamo and Non-Dynamo in Inertial-Range Turbulence 
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We investigate the Lagrangian mechanism of the kinematic "fluctuation" magnetic dynamo in 
turbulent plasma flow at small magnetic Prandtl numbers. The combined effect of turbulent ad- 
vection and plasma resistivity is to carry infinitely many field lines to each space point, with the 
resultant magnetic field at that point given by the average over all the individual line vectors. As 
a consequence of the roughness of the advecting velocity, this remains true even in the limit of zero 
resistivity. We show that the presence of dynamo effect requires sufficient angular correlation of 
the passive line- vectors that arrive simultaneously at the same space point. We demonstrate this in 
detail for the Kazantsev-Kraichnan model of kinematic dynamo with a Gaussian advecting velocity 
that is spatially rough and white-noise in time. In the regime where dynamo action fails, we also 
obtain the precise rate of decay of the magnetic energy. These exact results for the model are ob- 
tained by a generalization of the "slow-mode expansion" of Bernard, Gaw§dzki and Kupiainen to 
non-Ifermitian evolution. Much of our analysis applies also to magnetohydrodynamic turbulence. 

PACS numbers: 47.65.Md, 52.30.Cv, 52.35.Ra, 52.35.Vd, 95.30.Qd 



I. INTRODUCTION 



Turbulent magnetic dynamo effect is of great impor- 
tance in astrophysics and geophysics [l| . Many questions 
remain, however, about the basic mechanism of dynamo 
action, even for the kinematic stage when the seed mag- 
netic field is weak and does not react back on the advect- 
ing velocity field. Stretching of field lines by a chaotic 
flow is, of course, the ultimate source of growth of mag- 
netic field strength. Plasma resistivity rj in turn acts to 
damp the magnetic field. However, the dynamo cannot 
be understood as a simple competition between growth 
from stretching and dissipation from resistivity. For ex- 
ample, resistivity plays also a positive role in dynamo 
effect through the reconnection of complex, small-scale 
field-line structure [2]. 

In addition, random advection may not lead to field 
growth in the limit of vanishing resistivity. Consider, for 
instance, the kinematic dynamo model of Kazantsev Q 
and Kraichnan [J, with a Gaussian random velocity 
that is delta-correlated in time. In this model there is 
a dramatic dependence of dynamo effect on the spatial 
rugosity of the velocity, as measured by the scaling expo- 
nent < ^ < 2 of the spatial 2-point velocity correlation 
0. There exists a certain critical value ^, such that 
for ^ < kinematic dynamo effect exists only above 
a threshold value Pre of the magnetic Prandtl number 
Here k = ric/An is magnetic diffusivity 
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while ly is an effective viscosity associated to a "dissipa- 
tion length" £^ of the velocity, above which scaling holds 
with exponent ^ and below which the synthetic field be- 
comes perfectly smooth. In this regime of extreme rough- 
ness of the advecting velocity field, there is no kinematic 
dynamo even as k, ^ 0, if Pr < Pre- 

On the contrary, in the Kazantsev-Kraichnan (KK) 
model for smoother velocity fields with ^ > ^* there is 
a critical value Re„i^c of the magnetic Reynolds number 
Rem = UrmsL/n, whcrc Urms is the root-mean-square 
velocity and L is the integral length-scale of the fiuc- 
tuating velocity. See [1, 0], also In this smooth 
regime, small-scale kinematic dynamo leads to exponen- 
tial growth of the rms magnetic field, even for Pr — > as 
long as Rcm > Rem,c- The most natural correspondence 
of the KK model to the kinematic dynamo problem in 
real fluid turbulence is for the value ^ — 4/3, which is 
greater than the critical value = 1 in three space di- 
mensions. This correspondence would suggest that there 
is a critical magnetic Reynolds number for onset of kine- 
matic dynamo in actual fluid turbulence, but no lower 
bound on the magnetic Prandtl number. On the other 
hand, numerical studies of Schekochihin et al. [1, 
suggested that hydrodynamic turbulence is instead like 
the rough regime of the KK model and that the critical 
Prandtl number Pre for small-scale dynamo action tends 
to a finite, positive value as Rem — oo. Their latest in- 
vestigations now support the opposite conclusion, that 
the critical magnetic Reynolds number Rem,c tends to 
a finite, positive value as Pr o[lll,[T3- Considerable 
debate still continues, however, about the precise nature 
and universality of the observed small-scale dynamo. 

To resolve such subtle issues a better physical under- 
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standing is required of the mechanism of the turbulent 
kinematic dynamo. In our opinion, important ideas have 
been contributed recently by Celani et. al. 13]. They 
pointed out that the existence of dynamo effect in the 
KK model for space dimension d — 3 should be closely 
related to the angular correlation properties of material 
line- vectors. They considered the covariance at time t 
of two infinitesimal line-vectors that are advected start- 
ing a distance r apart at time 0. Celani et al. argued 
that this correlation vanishes as r decreases through the 
incrtial scaling range or as i ^ cxd, going to zero as a 
power {r/t^/^)'^ . Here 7 = 2 — ^ and Q = C(0 is the 
scaling exponent of a "homogeneous zero-mode" for the 
linear operator M.2 that evolves the pair correlations of 
line-elements forward in time. Celani et. al. further 
claimed that the transition between dynamo regimes in 
the KK model for d — i corresponds exactly to the value 
= 1 where C(?*) = [H. 

In this paper we shall further investigate these ques- 
tions. In the first place, we shall show that the claims of 
Celani et. al. [l^ are not quite correct. It will be shown 
here that the specific correlation function proposed by 
those authors docs not discriminate between dynamo and 
non-dynamo regimes. The scaling law which they pro- 
posed is valid, but holds over the entire range < ^ < 2 
with a different zero-mode and different scaling exponent 
than they had claimed. We shall show that a quite differ- 
ent correlation function of material line-elements is nec- 
essary to serve as an "order parameter" for kinematic 
dynamo. The crucial difference is that the quantity in- 
troduced here measures the angular correlation of ma- 
terial line-vectors that are advected to the same space 
point at time t. But still, why should there be any con- 
nection of roughness exponent ^ with dynamo action? 
Individual field lines ought to be stretched and their field 
strengths increased for all values of ^. We shall provide 
in this work a plausible physical explanation. Although 
individual lines may stretch due to chaotic advection, 
infinitely-many magnetic field lines will arrive at each 
point of the fluid due to diffusion by resistivity and the 
final magnetic field will be the average value that results 
from reconnection and "gluing" of field lines by resistiv- 
ity. We shall show that too little angular correlation leads 
to large cancellations in this resistive averaging, with the 
net magnetic field suffering decay despite the growth of 
individual field lines. 

We devote the remainder of our paper to a detailed 
study of the "failed dynamo regime" in the KK model 
for ^ < ^» and Pr < Pre- Part of our motivation is 
the speculation of that hydrodynamic turbulence 

at high magnetic Reynolds number but low Pr resem- 
bles this parameter range of the KK model. A bet- 
ter understanding of this regime may be useful to rule 
out its validity for hydrodynamic turbulence, based on 
astrophysical observations. We shall see, for example, 
that it implies a very rapid rate of decay of an initial 
seed magnetic field. Indeed, we show that in the KK 
model for rougher velocities the decay of the magnetic 



field is not resistively limited, with dissipation rate non- 
vanishing even in the zero-resistivity limit k — > 0, as long 
as Pr < Pre- There is a strong analogy with the anoma- 
lous decay of a turbulence-advected passive scalar, for 
which scalar dissipation is non- vanishing even in the limit 
of zero scalar diffusivitv(l3. [3]. The decay rate is in- 
stead determined by large-scale statistical conservation 
laws, associated to "slow modes" of the scalar evolution 
operator. We show here that the decay of the magnetic 
field in the rough regime of the KK model is determined 
in the limit k — > 0, Pr < Pre by the "slow modes" of 
the linear evolution operator for pairs of infinitesi- 
mal line-elements. We shall establish these results by a 
formal extension of the slow-mode expansion of Bernard 
et al. ^Td\ to the case of non-Hermitian evolution op- 
erators, which is presented in the Appendix. We shall 
furthermore determine all self-similar decay solutions of 
the magnetic field in the non-dynamo regime of the KK 
model, following for the passive scalar. Unlike the 
scalar case, however, determining the decay law of the 
magnetic energy requires an additional step of matching 
these self-similar solutions to explicit resistive-range so- 
lutions. We shall use these results to discuss the physical 
mechanism of kinematic dynamo, and, in particular, to 
relate our dynamo "order parameter" to the process of 
"induction" by a spatially uniform initial magnetic field. 
As we shall see, considerable insight can be obtained into 
the inner workings of the small-scale dynamo by consid- 
ering also the situations where it fails. 

II. THE KRAICHNAN-KAZANTSEV DYNAMO 
AND CORRELATIONS OF LINE-ELEMENTS 

A. The Kinematic Dynamo 

The evolution of the passive magnetic field B(x, i) is 
governed by the induction equation 

9tB + (u.V)B- (B.V)u = kAB, (1) 

where u = u(x, t) is the advecting velocity field and k is 
the magnetic diffusivity. The magnetic field is taken to be 
solenoidal, assuming there are no magnetic monopoles: 

V-B = 0. (2) 

Notice that this condition is preserved by the evolution 
equation ([T]) if it is imposed at the initial time to — 0. 
We have also assumed above that the advecting fluid is 
incompressible so that 

V-u = 0. (3) 

For simplicity, we shall only consider this case hereafter. 

For an incompressible fluid, one can represent the solu- 
tion of the induction equation by a stochastic Lagrangian 
representation of the following form: 

B(x, t) = E [Bo(a).Vax(a, t) |a=a(x,t)] • (4) 
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See [Tj, [T3| ■ Here a(x, t) are the "back-to-label maps" 
for stochastic forward flows x(a, t) solving the SDE 



dx(a, t) = u(x(a, t),t)dt + V2KdW{t). 



(5) 



E denotes average over the realizations of the Brownian 
motion W(t) in Eq. ([5]). For the ith component of the 
magnetic field we can write 



B*(x,i) = E 



5o(a)^(a' i)la=a(x,t) 



(6) 



= Jd''aB^{a)Fl{a,0\K,t;u) 
where we have defined 



f^,!(a,0|x,t;u) = E 



^MS^' (a-a(x,i)) 



(7) 



These results may be used to represent the equal-time, 
2-point correlation of the magnetic field, averaged over 
the stochastic velocity u and the random initial magnetic 
field Bo : 

(B*(x,i)BJ(x',<)) = /d''a/dV(B^-(a)S^(a')) 

Ffc',(a,a',0|x,x',i). 

We have assumed that u and Bq are statistically inde- 
pendent and we also have defined 

Flf,(a,a',0|x,x',t) = (F^,(a, 0|x, t)/'/ (a', 0|x', t)\ . (8) 



For statistically homogeneous velocity and initial condi- 
tions, with C'^r,t) = (B*(x,i)BJ (x',t)) for r = x - x', 
we obtain 



(9) 



with 



EE' 



^(a + p,t)^(a,i) 



(a(x + r,f) -a'(x, t) - p) 



a— a'(x,^) 

(10) 



Here the prime / denotes a second Brownian motion 
W'(i) statistically independent of W(i). Equation ^ 
was introduced by Celani et al. [l3| and heavily exploited 
in their analysis of the magnetic correlation. 

Another closely related propagator was introduced by 
related to infinitesimal material line-elements, which 
evolve according to the Lagrangian equation of motion: 

DM = (M-V)u. 

Note that the positions of the line-elements are assumed 
to move stochastically according to (O, so these are not 



quite "material lines" in the traditional sense when k > 0. 
The exact solution of the above equation for i > is 

bi{t) ^ «(0)-Vax(a,i), 

with x(a, i) solving ([5]). Taking initial line-elements 

'^^fc(O) = ^fc' '^^^■'(0) = ^\ starting at positions a, a' dis- 
placed by r = a' — a, one may follow to define for 
statistically homogeneous turbulence 



x'(i)-p)) 



= EE' 



9x* , , dx''' 



(a,t) 



X(5'^(ic(a + r,t)-x'(a,i)-p)^ (11) 

If we make the change of variables a i—s- x in the argument 
of the delta function of eq. fTUI) . then the jacobian of this 
transformation of variables is 1 due to incompressibility. 
Therefore, one finds by comparison with (jlip that 



FZ(.P,0\r,t) = EE' 



Ox^ 



(a,i) 



(i(a + p, <) - i'(a, t) - r) 



equating the two propagators under interchange of argu- 
ments. 

In our work below an important role will also be played 
by the covariant vector given by the gradient G = V0 of 
a passive scalar 9. The gradient satisfies the equation 



dtG -f (u-V)G + (Vu)G = kAG. 



(12) 



which is dual to the equation ([T]) for the contravariant 
vector B [l8|. The above equation preserves the con- 
dition G — V6' if this is imposed at time = 0. A 
stochastic Lagrangian representation also exists for the 
solution of this equation. Solved forward in time with 
K > this representation involves the matrix Vxa(x, t). 
However, taking k, —k in (|12p and solving backward 
from time t > to time yields the representation for 
the i component of the gradient field: 



Gz(a,0) 



dx'' 



(a, i)Gfc(i(a,t),i) 



= J d'^x Gk{x, t)F^ (a, 0|x, i; u) 



For statistically homogeneous velocity and initial condi- 
tions we introduce the 2-point correlation of the gradient 
field, g.,j{p,t) = (G,(a,0)Gj(a',0)) with p = a - a'. By 
the same arguments as previously 
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gy(p,0) = d''rgM{v,t)F^Ap,Q\v,t) 
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= J d^rgki{r,t)Fl^/{r,t\p,0), (13) 

for positive times t > Q. 

We shall generally avoid using the geometric language 
of differential forms and Lie derivatives in this paper, but 
a few brief remarks may be useful. For those unfamiliar 
with this formalism, a good introductory reference is p^ . 
The magnetic field B discussed above is a 1-form, which 
is more properly represented by its Hodge dual B*, a 
{d — l)-form. The equation ([T]) for k = is equivalent to 
(9tB* + £B* = 0, where C is the Lie derivative. The Lie 
derivative theorem thus implies that ([T]) for k = satisfies 
an analogue of the Alfven theorem, with conserved flux 
of B through (d— l)-dimensional material hypersurfaces. 
On the other hand, the field G is a proper 1-form and 
equation l|12p for k = is equivalent to 5(G + CG ~ 0. 
The Lie derivative theorem thus implies that integrals of 
G along material lines are conserved for k = 0. Either of 
these equations could be regarded as a valid generaliza- 
tion of the kinematic dynamo problem to general space 
dimension d. The non-gradient solutions G of (fT^ are 
generalizations of the 3-dimensional vector potential and 
the "magnetic flux" is represented by their line-integrals 
around closed material loops. For d — S the non-gradient 
solutions of are in one-to-one correspondence (up to 
gauge transformations) with the solenoidal solutions of 
([T]), by the familiar relation B — VxG. 

We note finally that all of the results in this section 
hold for any random velocity field that is divergence- 
free and statistically homogeneous. Thus, they apply 
not only to the Kazantsev-Kraichnan model discussed in 
the following sections, but also to the kinematic dynamo 
problem for hydrodynamic turbulence governed by the 
incompressible Navier-Stokes equation. 

B. White-Noise Velocity Ensemble 

In the Kazantsev-Kraichnan (KK) model [1, 0, [E] the 
advecting velocity u(x, t) is taken to be a Gaussian ran- 
dom field with zero mean and second-order correlation 
delta function in time, given explicitly by 

(u'(x, i)u^ (x', t')) = [DoS^^ - S''{r)] S{t - t') (14) 

with r = x' — X. Under the incompressibility constraint 
diS'^^{r) = and supposing S^^{r) scales as ~ for 
<C r <C i, one deduces for that range that 

S'^{r) = Dir^ [i^ + d- 1) S'^ - ^f'f^ (15) 

where f* = r'^/r. Define "viscosity" v ~ Di^t- Be- 
low we shall consider especially the limit v, k ^ with 
v < PtcH. This is the non-dynamo regime in the limit of 
infinite kinetic and magnetic Reynolds numbers. One of 
our main objectives is to understand better the geomet- 
ric and statistical properties of this regime which lead to 
the failure of small-scale dynamo action. 



In addition to the properties of statistical homogeneity, 
stationarity and incompressibility, the white-noise veloc- 
ity ensemble is time-reflection symmetric. This implies 

F^,(r,-t|p,0) = F;-',(r,t|p,0) 

and the similar property for i^^;^ (p, t|r, 0) . Combined 
with the other symmetries, this imples that 

= F%{r,-t\p,0) 

= FZ{r,t\p,0). 

The first line follows from incompressibility, the second 
line is due to time-translation invariance and the last 
equality follows from time-reflection symmetry. Thus the 
two propagators are adjoints in the KK model. 

Time-reflection symmetry has also an important impli- 
cation for the evolution of the gradient field correlation. 
Note that time-translation invariance implies that equa- 
tion can be written for t > as 

g^,{p,~t)^ J d'irgM{r,0)Ftf{r,0\p,-t). 

Then time-reflection symmetry implies further that 

g^J{p,t)^ J d^rgki{r,0)F^/{r,0\p,t) (16) 

for t > 0. Compare with equation ^ for the magnetic 
correlation. We see that the i^-propagator in the KK 
model evolves forward in time the gradient correlation. 

The most important property of the white-noise model 
is its Markovian character, which implies that time- 
evolution of correlations is governed by 2nd-order dif- 
ferential (diffusion) equations. E.g. the n-th order equal 
time correlation function = (Ha^i (xa, i)) 

satisfies an equation of the form dtCn = M.nCn- Expres- 
sions for the general n-body diffusion operators A^„ can 
be found i nj20i] , which can be obtained using Ito formula 
as in Ref. [7[ or, equivalently, by Gaussian integration by 
parts. In the limit z^, k ^ all of these operators for gen- 
eral n become degenerate (singular) and homogeneous of 
degree —7 with 7 — 2 — ^. Below we shall mainly consider 
n ~ 2 and thus write simply M. for M.2- However, many 
of our considerations carry over to general n, as will be 
noted explicitly below. Following the notations of p^ . 
we write for n — 2: 

a,r^(r,t) = [M(r)];^C'"'(r,i) (17) 

with 

(18) 

The notation A^(r) indicates that di = d/dr"^. Note 
that equation ^\7\ has an invariant subspace satisfying 
diC-i = djC'^ = 0. 
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It follows from (|T7|) and ([9]) that the F-propagator is 
the heat kernel of the adjoint operator 

(19) 



satisfying 



dtF^,{p,0\v,t) = [M*ip)]ll_F^^{p,0\r,t) 

= [M{v)]lFZ{pAr,t). (20) 

The propagator F is thus the heat kernel of Ai: 

9tJ^«(r,0|p,<) = [Mirt;^Frf{r,0\p,t) 

= [M*ip)r^jF^^{v,0\p,t) (21) 

Because of the homogeneity of the operators Ai and A4* 
in the i/. k — > limit, satisfies the scaling relation 

J^f (Ar, 0|Ap, X-'t) = X-'^Ftf (r, 0|p, <) , (22) 



with an identical relation for the _F-propagator. 

Finally, it follows from (fT6|) that the gradient correla- 
tion satisfies 



dtg^,{p,t)^[M*ip)r,JgpMt)- 



(23) 



This equation has an invariant subspace of solutions of 
the form Qij = —didjQ for a scalar correlation function 
0(r, i). Celani et al. [l^ have also introduced the quan- 
tity 



QM(r,t) = / d^pFl\[p,t\r,Q) = {5lk{t).5l',{t)) , (24) 



where on the righthand side the line-elements are ini- 
tially unit vectors ^£fc(0) = e^, (5£^(0) = starting at 
positions displaced by r. This quantity measures the an- 
gular correlation of the material line-elements at times 
t > 0, as well as their growth in length. It follows from 
([21) that this quantity in the KK model satisfies 



(25) 



with initial condition Qfc^ (r, 0) = Ski, as already noted in 
[T3| . This equation is identical to ([^5]) for the gradient 
correlation and, furthermore, Qfc£(r, 0) = — 9fc9f9(r,0) 
with e(r,0) = -(l/2)r2. Thus, Q is of gradient type. 

In this work we restrict ourselves to conditions of sta- 
tistical homogeneity, isotropy and parity invariance for 
all stochastic quantities. Thus, we can write the 2-point 
correlation function of the magnetic field as 

^ Cl (r, tyr^ + Cn (r, t) [5'^ - f'f^ ) (26) 

where f* = r'^ /r. Cl and Cn are the longitudinal and 
transverse correlations, respectively. With the form of 
the velocity correlation in psp . the evolution equation 
(IT71) reduces to two coupled equations for Cl and Cn- A 



lengthy but straightforward calculation gives 
dtCL 

= Dir^ {{d-l) drrCL + (d + 1) (d - e - 1) ^drCL 

+ id-l) [e-i-2{d-l)] ^Cl 
+ {d-l)[{d+l)C + 2{d-l)]^CN}. 

and 
dtCN 



(27) 



'^-drC 



^ Dir^ i^{d^l)drrCN + [(d+l)^+(d-l)- 
+ [{d+l)e + {d^-5)(~2{d~l)] ^Cn 
+ (e-2) {^-l){d + ^-l)^CL} , 

(28) 

respectively, when VjK 0. For solenoidal solutions, 
such as for the magnetic field, it is easy to show that 
the two correlations are related by 



N 



Cn = Cl 



1 



d-l 



rdrCL- 



(29) 



For example, see [1^ |2l|]. The solutions satisfying this 
relation form an invariant subspace, with the evolution 
reducing to a single equation for Cl'- 

dtCL = Dir^ [id - 1) drrCL + (2? + d^ ~ l) ^drCL 

+ {d-md + OT^CL] - 

(30) 



In the same manner, the general solution of ([23]) or 
(1^ may be decomposed as 

Gi] = GL{r,t)fifj + GN{r,t){5i-j - hfj) (31) 

under assumptions of homogeneity, isotropy and 
reflection-symmetry. Then Gn and Gl satisfy the fol- 
lowing coupled equations for i^, k 

dtCL 



= (d - 1) Diri {drrCL + (3^ + - 1) ^drCL 

+ (3e + 2d-2) {^-1)^Gl 

+ id + ^-l)i^-l){^~2)^GN} 



(32) 



'N 



(d-3)e + (d-i)' 



and 

dtCN 

= Dir^ ^^{d-l)drrG 

+ [(d+l)^ + 2{d-l)]^GL 
+ [{d- 1)^-2] {^ + d-l)^GN}, 

respectively. Gradient solutions satisfy the constraint 

Cl^Gn + rdrCN (34) 



N 



(33) 
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where Gn = —^drQ in terms of the scalar correlation 
function O. In this invariant subspace of solutions the 
dynamics reduces to a single equation for Gn- 

dtGN = id-l)Diri[drrGN + {2^ + d+l)^drGN 

(35) 

C. Line- Vector Correlations 

Kinematic dynamo effect is due ultimately to the 
stretching of magnetic field lines as they are passively 
advected by a chaotic velocity field. However, the prop- 
erties of infinitesimal material line-elements in the KK 
model are, at first sight, counterintuitive in this respect. 
In order to discuss stretching of individual lines, it must 
be understood that the velocity field is smoothed at very 
small scales < £i^. The inertial-range velocity structure 
function in (jlSp crosses over to a viscous-range form 

S'y (r) = DiEi'^r^ [(d+l) 5'^ - 2r'f^] (36) 

for r ^ (-u- The growth of line-elements in such a smooth 
velocity field, white-noise in time was derived by Kraich- 
nan |22l | to be exponential 

{5f{t)) « e^^* (37) 

with the Lyapunov exponent 




= Diii-'^d{d-l). (38) 

See also [H [13, HI]. The "material" line-elements of 
relevance to the kinematic dynamo are subject to an ad- 
ditional Brownian motion proportional to \/2k in ((5|). 
However, this effect of molecular diffusivity k corresponds 
just to changing the constant Z?o in the Kraichnan-model 
velocity covariance (|14p to i^o + 2k. Since only velocity- 
gradients enter (|38p . this result for the Lyapunov expo- 
nent still holds in the presence of molecular diffusivity. 

It follows from p8|) that line-stretching is greater for 
smaller ^ and smaller v. This may seem a bit perplexing, 
because the dynamo fails for ^ too small, in the range 
< ^ < In that regime, there is no dynamo ac- 
tion for V < PrcK, despite the fact that the stretching 
rate becomes larger as v decreases. The turbulent kine- 
matic dynamo cannot be understood as a simple "com- 
petition" between stretching and diffusion. What, then, 
can account for the presence of the dynamo in the range 
^* < C < 2 of smoother velocities in the KK model and 
failure of the dynamo in the range < ^ < ^, of very 
rough velocities? An intriguing suggestion has been made 
by Celani et al. [l3] that the existence of dynamo effect 
can be characterized by the angular correlations of mate- 
rial line-elements. They proposed the function Q defined 



in p4|) as an "order parameter" for the dynamo transi- 
tion. As we shall demonstrate below, the principal con- 
clusions of [l3] about Q are erroneous and this quantity 
does not discriminate between dynamo and non-dynamo 
regimes in the KK model. However, our discussion will 
lead us to identify a different correlation property of in- 
finitesimal material line-vectors, which can indeed serve 
as an "order parameter" for the dynamo. 

The principal claims of [13] were as follows. First, in 
the non-dynamo regime for < ^ < ^* with rough veloc- 
ity, Q exhibits a scaling law of correlations: 

Qfc,(r,0^ (const.) Zki{r), (39) 

for r ^ inm{{Dity/'^ , L}. For finite z^, k, this re- 
lation holds in the inertial-convective range of scales 
max{£„,£i,} <^ r ^ L, with = {n/DiY^^ (assuming 
that V < PrcK). This is an example of "zero-mode domi- 
nance" Thus, the quantity Zkiir) is a homogeneous 
zero mode of the operator M* ^ satisfying [A^*]^'Zpg = 0, 
with exponent C(0 > for < ^ < Intriguingly, it 
was found that ^ = ^ + 2, where C is the scaling exponent 
of the zero mode of M. which was shown in 26] to dom- 
inate in the magnetic correlation of the KK model for 
the same parameter range. It was furthermore claimed 
in [13] that C(C*) = 0- For 2 > ^ > on the contrary, it 
was argued that M* develops point spectrum and that 

Qm{v, t) ^ (const.)e^«*£fc,(r), (40) 

where Eq is the largest positive eigenvalue of Al* and 
£k(.{^) is the corresponding eigenfunction. Eq is numer- 
ically equal to the dynamo growth rate. Notice in the 
limit Pr < 1 that Eo (xl/t^ = {Df/K'*)^^^ [1], so that 
K must be kept nonzero (but with Rem > Rem,c)- Thus, 
the "material lines- vectors" in Q are advected by velocity 
u subject to Brownian noise proportional to y/K. The ap- 
propriate terms proportional to k must then be included 
in the diffusion operators M and M* jl3| . 

We shall show that the zero-mode dominance rela- 
tion (|39p does hold for Q but with a different zero- 
mode and different scaling exponent C than that claimed 
by [l3|. Furthermore, the scaling relation (I39|) holds 
over the whole range < ^ < 2, assuming only that 
max{^„,£^} <c; r <C L, with an exponent C(C) = ~C < 
which does not exhibit any qualitative change at the dy- 
namo transition = ^* . The exponential growth relation 
(UOj) does not hold for the quantity Q anywhere over the 
range < ^ < 2, even if Pr > Pre and Rcm > Rem^c- 

D. Zero-Mode Analysis 

The basic tool of our investigation is a generalization 
of the slow-mode expansion of Bernard et al. [iGj . Those 
authors derived such an expansion for the propagator or 
heat kernel P(r, t\ro, to) that describes the evolution of a 
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passive scalar in the Kraichnan white-noise velocity en- 
semble with covariance (|14p . However, the derivation of 
[l^ was, in fact, axiomatic and applicable to the prop- 
agator for any non-positive, self-adjoint operator, with 
absolutely continuous spectrum and homogeneous of de- 
gree —7. This derivation showed that the slow mode ex- 
pansion follows from assumed meromorphic properties of 
the Mellin transform of the propagator and Green's func- 
tion of the operator. In the Appendix of this paper, we 
generalize this axiomatic derivation to the case of the 
non-Hermitian operators M and M* , in the non-dynamo 
regime where both have absolutely continuous spectrum. 
We refer the reader to the appendix for details and here 
just state the essential results. 

The operators Ai and A4* have two types of homo- 
geneous zero-modes, regular and singular. The regular 

zero modes are denoted Z^^) and Z^"'' for a — 1,2,3, 
respectively, with exponents Ca and Ca whose real parts 
increase with a. These are ordinary functions which sat- 
isfy the conditions MZi^^) = and M*z'"°'^ = globally. 

The singular zero modes W(a) and W^""^ for a = 1, 2, 3, 
instead have exponents uja and ZJa whose real parts de- 
crease with a for a — 1,2,3,..., respectively. These are 
distributions which satisfy the conditions A4W(^a) — 

and A^*W*'°'' = only up to contact terms. The scaling 
exponents of the two sets of zero modes are related by 



(41) 



Above each regular zero mode lies an ascending tower 

of slow modes Z(^a,p) and z^""'^^ homogeneous of degree 

Ca,p = Ca +7P and Ca,p = Ca+lP: respectively, satisfying 

A4Z(^a,p) = ^(a,p-i) and A^*z' ''''' = Z^ "* for p — 

1,2,3,... with Z((j o) — ^(a) and Z*' — z'' \ In addi- 
tion, below each singular zero mode is a descending tunnel 

of self-similar decay solutions M^(a.p) (r, t) and W^"^'^"^ (r, t) 

satisfying dtW(^a,p) = -^T^Ca.p) and dtW^""'^^ = 

M*W^"'^^ with T4^(,^p)(Ar,A^t) - A-»-(P+i)W(,,p)(r,t) 

and W^''"''\Xr,X''t) = X'^-^-^P+^'^''W^''"''\r,t). These are 
related to the singular zero modes by MW(^a.p-i) = 

W(^a,p} and M*W^''''"^^ = "W^^"'^^ for p = 0, 1,2, ... with 

Ty(„,_i)(r,0) =M/(„)(r) and W^''^''\r,0) =W^''\r). 

In terms of these quantities there are short-distance 
expansions, for A ^ 1, both for the F-propagator 

rT77(<^.P) / 



F^,iXr,t\p,0) ^ ^ X^-^^-^PZli^rW^rip^t)]*, 

a,p>0 



and for the _F-propagator 



(42) 



F;',(Ar,t|p,0) 



a,p>0 



XC.+^PziY\r)[WlJp,t)] 



See the Appendix for the details of the derivation. Note 
that these asymptotic series are generally dominated by 
their leading terms for a — 1 and p = 0, corresponding 
to the regular zero mode with scaling exponent of small- 
est real part. Of course, the leading term may give a 
zero contribution for various reasons and then sublead- 
ing terms will dominate instead. In order to make use of 
this expansion we must calculate explicitly the homoge- 
neous zero modes of M and M* . 

To find the isotropic and scale-invariant zero-modes of 
A4 we substitute into ([27| and ((28|) the forms 

CL^A^r'^, CN = ANr^ 

giving the matrix equation 



Mll 


Mln ' 




■ Al ' 




■ 0" 


Mnl 


Mnn _ 




_An _ 








with 

Mll = {d-l)a{a-l) + {d-l){d-^-l)a 
+ (d-l)(e'-e-2rf + 2) 

MLN^{d-md+l)i + 2{d-l)] 

MML^{^-2){^-l){d + ^-l) 

Mnn = {d-l)(T{(j-l) + [{d+l)£, + {d~lf]a 
+ {d+l)^'^ + {d^ -b)^-2d+2. 

Calculating the determinant 



Mll Mln 
Mnl Mnn 
[(d-l)a2 + 



= (d-l)(CT-2)(o- + d-2)x 

{d^-d + 20a + {d-l)i{d + C)\, 
one finds that the scaling exponents a are 

1/2 



Ci 



UJ2 



d__ 

' 2~ d 



1 



d 

2 



d-l 



1-4^ 



(d-2)(d + e- 1) 



1-4^ 



d{d - 1) 



(d-2)(d-t-e- 1) 



1/2 



rf(d- 1)2 

C2 = 2, uji = 2-d. 

Note that the set Ci i ^2 correspond to the invariant sub- 
space of solenoidal solutions, as may be verified by sub- 
stituting the scaling ansatz for Cl into ([30]) . The expo- 
nent Ci coincides with that found by Vergassola [1^ to 
dominate in the magnetic 2-point correlation for a forced 
steady-state at high magnetic Reynolds number and zero 
Prandtl number. 

To find the isotropic and scale-invariant zero-modes of 
A^* we likewise substitute into (f32|) and ([33|) the forms 



(43) 



Gl 



Ajr° 



G 



N 



AnV 
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giving the matrix equation 



Mll Mln 





Al 

















with 



+ (d-l)(3e + 2d-2)(e-l) 
Mln = {d-l){d + ^- - m - 2) 
Mnl = (d+l)e + 2(d-l) 



Mnn = (d- IM^- 1) + + (rf- 1)2)ct 

+ ((d-l)e-2)(d + e-l). 

Calculating the determinant 
Mll Mln 



M M -{d-l){a + S){a + i + d)y. 

M NL iVl NN 

[{d - 1)^2 + ((d - 4)((i - 1) + 2^{d - 2))W 

+2(d-2)(d + e-l)(^-l)], 



one finds that the scaling exponents a are 



A-d 2-d d 



1-4^ 



(d-2)(d + g- 1) 
did -If 



1/2 



A-d 



2-d d 
2~ ^ d-r^2 



1-4^ 



(d-2)(d + ^- 1) 



-,1/2 



d{d-lf 



Note that the set C,i , cJ2 corresponds to the invariant sub- 
space of gradient solutions, as may be verified by substi- 
tuting the scaling ansatz for Gat into ([55]) . These expo- 
nents are directly related to those for the passive scalar 
discussed in [l^. In the isotropic sector of the scalar, the 
regular zero mode has exponent (constants) and the sin- 
gular zero mode has exponent j — d. The first vanishes 
after taking gradients and is replaced by its lowest-lying 
slow mode, with exponent 7. Taking two space deriva- 
tives reduces these exponents by 2, yielding 7 — 2 = — ^ 
and (7 — d) — 2 = —{d + ^). The exponents C2;'^i coin- 
cide (for d = 3) with the set which were claimed in [l^ 
to determine the scaling of Q. 
It is easy to check the relations 

'^a + Ca = -d + l, UJa+Ca = -d + J 

for a — 1,2, consistent with the general result (|¥T|) . We 

see that, at least in the isotropic sector, the solenoidal 

scaling solution W^(2)(r,t) is associated in the slow-mode 

(2) 

expansion to the non-gradient zero-mode Z (r) and the 



(2) 

gradient scahng solution W (r, t) is associated to the 
non-solenoidal zero- mode 2'(2)(r). In fact, this is true in 
general, as we now show. Take any solenoidal scaling 
solution W(^a.p){p, t)- Then it follows from the propagator 
relation © and the scaling property of the singular slow- 
mode W^(a,p) (r, 0) that 



d'^r 



W; 



kl 



(r,0) F^,(Ar,i|p,0) = A^"+^^' 



This can only be consistent with the slow-mode expan- 
sion dlSl) of F for A < 1, if 



W^f4)(r,0) ZlY'(r)^0 



Since M^(a,p) is solenoidal, then Z ' must be non- 
gradient. Otherwise the integral will vanish, because the 
solenoidal and gradient subspaces are orthogonal. An 
identical argument shows likewise that any gradient scal- 
ing solution is associated in the slow-mode expan- 
sion of to a non-solenoidal zero-mode Zf^^^^py 
As should now be clear, however, the relation 



can- 



<2) 



. Since Q is of gradient 

(1) 



not hold with C = C2 ^-^^ Z = Z 

type, its evolution is described by (I35|) which has as 
its only regular zero mode with scaling exponent Qi. We 
shall now verify this directly from the definition (|24|) of 
Q, by means of the slow-mode expansion. We use first 

the adjoint relation F^^(p,i|r,0) = ^^^^(r, t|p, 0) and the 



homogeneity relation ((22|) for F to write 

Q«(r,t) = j d''pFZiXr,l\p,0) 

with A = r/{Dit)'^/-' and p = p/iDitfl^. Then using 
dU) for A < 1 gives 



a,p>0 



Notice, however, that the space integral vanishes for all 
Wra,p) in the solenoidal sector (e.g. see This follows 



for any solenoidal correlation C , from 

C"(r,t) = dkdeAke{r,t) - AAkk{r,t), 



(44) 



where Ake is the correlation of the vector potential A. 
Note that, in general dimension d, B is a 1-form and 
A is a 2-form, related by the codifferential B = (5A[4oI|. 
Since the solenoidal solutions W^(a,p) are associated in the 

expansion to the slow modes Z^"^'^^ outside the gradient 
sector, all of these terms drop out in Q. The result is the 
same as the slow-mode expansion carried out entirely in 
the gradient sector, with the leading term 



{ki 



(r,t)~C2 



(45) 
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for C2 = / d'^pW^2 o)iPj^)- This is the correct relation 
replacing relation ([59)1 claimed in [l^l ■ 

This same relation may be verified by appealing to the 
results of Eyink and Xin 14] on the self-similar decay of 
the passive scalar. Those authors found that there is a 
universal form of the self-similar decay solutions for the 
scalar correlation function at short distances: 



x{t) 



(46) 



for r < (Dii)!/'^. See equation (3.21). Here 

x{t) — — {1/2) {d/dt)'d'^(t) is the scalar dissipation rate. 
For general initial data with power-law decay of correla- 
tions in space, 8(r, 0) ^ r~", the decay rate is given by 
-d^it) - ^-"/^ at long times [H. Since Sfe^(r,0) = Sm 
corresponds to Q{r,0) — — (l/2)r^ with a = —2, we re- 
cover from Qfcf(r, t) = —dkdiQ{r,t) and eq. ip^l) exactly 

the relation (|45)l . with z'"^^ (r) = 5m — ^rkrg. The latter 



result may be verified from equation p4p by substituting 
Gn — Ant^^- This alternative derivation of (j45ll makes 
clear its validity over the whole range < ^ < 2 and not 
just < ^ < There can be no exponential growth 
relation for Q, such as relation (^0]) proposed in 13]. It 
is true that the operator A4* must have a positive eigen- 
value whenever Ai does so. However, the corresponding 
eigenf unctions must lie in the non-gradient sector. An 
exponential growth for Qke{j^,t) = —dkde'd{r,t) would 
require an exponential growth for the scalar correlation 
function <d{r,t), which does not occur. 



E. A Dynamo Order-Parameter 

Based on the previous discussion, we now will propose 
an alternative definition of a line-correlation which can 
serve as an "order parameter" for the dynamo transition. 
Clearly, one should not integrate Fll{p,t\r,0) over p, as 
this eliminates the solenoidal sector. We propose instead 
to set p = 0, defining: 

7^fe,(r,^) = Fll{O,t\v,0) 

= (%(t).(54(i)5''(x(t)-x'(t))), (47) 

where, as before, the two line elements are started with 
(54(0) = efc, (54(0) = ee, and x'(0) - x(0) = r. Because 
of the delta- function, TZke{r,t) measures the growth in 
magnitude and the angular correlation between those ma- 
terial line- vectors which arrive, stretched and rotated, at 
the same point at time t. Just like the quantity Q defined 
in [1^, TZ satisfies also the equation 

at7^fc,(r,^) = [M*{r)]ljnpg{r,t). 

However, it has the initial value 

7^H(r,0) =(5M(5'^(r), 

which is non-gradient, unlike for Q. Thus, TZ should expe- 
rience exponential growth like (|40p in the dynamo regime 
for 2 > ^ > . 



The time-dependence in the non-dynamo regime for 
< ^ < ^* can be obtained from the slow-mode expan- 
sion of F'J^(Af, 110,0) with A = r/{Dity^^. One obtains 
for A < 1 that 



a,p>0 



Thus, TZ exhibits a power-law decay in time, with the 
dominant terms given by the two isotropic zero modes 



7^fef(r,^) 



±hii^(l) 



CiiDit) 



M V- 



Zui^) (48) 



In all dimensions d the exponent C2 > for < ^ < 1, 
becoming negative for ^ > 1. Thus, the first term with 
— — ^ dominates for lower ^ values. There is a critical 
dimension d^ = 4.659, given by the positive real root of 
the cubic polynomial d^ — 8d^ + 19d— 16, above which it 
instead true that (^2 < Ci when ^ > £,c with 



v/(d2 -3d + 4)2 + Sjd - l)2(d - 2) - (d^ -3d + 4) 
2(d- 1) ■ 




Note that S.c < where the dynamo transition occurs. 
The latter value 0] 



= (rf - 1 



is the point at which C2 develops an imaginary part and 
the slow-mode expansion above breaks down. 

For exponents < ^ < 2 in all the integer dimensions 
2 < d < 9, the power-law decay is replaced by exponen- 
tial growth 



(r). 



(49) 



proportional to the eigenfunction £fc£(r) of A^* with 
largest eigenvalue Eq. To demonstrate this, it is enough 
to show that the initial condition TZk£{r,0) gets a non- 
zero contribution from the eigenfunction £ke{r). We may 
represent this initial state by an expansion 



7^fef(r,0) 



C'd^/cf (r), 



where £fc^(r) is the eigenfunction of A^* with eigenvalue 
Ea- Note that for the continuous spectrum, this is a gen- 
eralized eigenfunction expansion where the sum over a 
is a continuous integral and S/.^ are distributions, not 
square-integrable functions. The expansion coefficients 
are given by 



a = {£c„TZm - J d''r£l\v)TZkt 



(r,0) 



where are the eigenfunctions of M with the same 
eigenvalue These form a bi-orthogonal set with the 
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eigenf unctions E^.^ of A^* [4l|. The coefficient Co cor- 
responding to the eigenfunction £f.^ ~ Ekt is non-zero 
because 

where £'^'^(0) 7^ is (twice) the energy in the normalized 
dynamo state. 

Thus, unhke the quantity Q proposed in [l3|, the hne- 
correlation TZ defined in (|T7)l satisfies the exponential 
growth relation ()49p in the dynamo regime and power- 
law scaling (|48p in the non-dynamo regime. It would 
be of great interest to determine the spatial structure 
of the eigenfunction ^^(r). Of course, this function 
must be of non-gradient type. It is known Q that 
the trace of the dual eigenfunction £{r) = £"(r) ex- 
hibits stretched-exponential decay of the form £{r) oc 
— exp (— /^Cj^/^k)^^^) for r ^ £f^. A similar behaviom for 
£{r) = £kk{j^) can be checked to be consistent with the 
dynamical equations, but a more careful investigation is 
required. This will be pursued elsewhere. 

A quantity with even simpler geometric significance 
which might also serve as an "order parameter" is 



7^(^) 



d''rnkk{r,t) = {S£{t)-5£'{t))o. (50) 



This is the covariance of two material line-elements which 
started at any relative positions as identical unit vectors 
at time and which ended at the same point at time 
t. The notation denotes the conditional expectation 
over material lines which end at zero separation. Clearly 
TZ{0) = 1. However, its time-dependence is undetermined 
by our present considerations both in the dynamo and in 
the non-dynamo regimes. We cannot argue that TZ(t) 
decays as a power in the non-dynamo regime, because 
the slow-mode expansion applies only for r ^ (Dity^'' 
whereas the definition of Ti-it) involves an integral over 
all r. We also cannot conclude that TZ{t) grows expo- 
nentially in the dynamo regime, because this requires 
the condition/ d'^rSj-iir) ^ 0, which needs to be shown. 
However, we shall see in the next section that TZ(t) has 
a direct interpretation in terms of the turbulent decay of 
an initially uniform magnetic field and we shall determine 
its time-dependence in the non-dynamo regime. 



III. DECAY OF THE MAGNETIC FIELD 

We now consider in detail the problem of the turbulent 
decay of the magnetic energy {B'^{t)) in the non-dynamo 
regime of the KK model, for < ^ < ^, and Pr < Pre- 

A. Discussion of the Convective-Range Decay Law 

We begin by giving a simple, heuristic explanation of 
the decay law of the magnetic field, for generic initial data 



of the magnetic field with rapid decay of correlations in 
space. The fundamental observation is that the zero- 
modes a = 1,2,3,... of the adjoint operator M* 
give rise to statistical conservation laws in the evolution 
of the 2-point correlation. 



Ja{t) 



d''rZ[f{r)C\v,t), 



which satisfy 



{d/dt)Ja{t) = J d''rZ[fir)-[Miv)]%C'"^{r,t) 

= J d'r [M*{r)]%z[f{r) ■ C^"'(r, t) = 0. 

Note, however, that only the non-gradient zero modes 
lead to non-trivial conservation laws, because of the or- 
thogonality of solenoidal and gradient correlations. The 

(2) 

leading-order zero-mode is thus Z found in the previ- 
ous section, namely. 



with 



A 



{2) 



(e-2)[(d-l)C2 + (rf-3)(d + e-l) 



= (rf+l)C + 2(d-l). 

This zero-mode coincides with that found by Celani et 
al. |13|] for d — 3. 

The corresponding conserved quantity J2(i) plays the 
same role in the turbulent decay of the magnetic field 
as played by the "Corrsin invariant" in the decay of the 
passive scalar [3; 13 ■ Assume, in fact, a self-similar 
decay law for the magnetic correlation 

C^r,t) = lt^{t)r^r/L{t)). 

The length L{t) is a large-distance correlation length or 
"integral length" of the magnetic field. The quantity ti (t) 
is a measure of the magnitude of the magnetic fluctua- 
tions at scale L(t), which we term the magnetic ampli- 
tude. Just as for the scalar, the growth of the magnetic 
length-scale L{t) can be obtained dimensionally from 



1 d 
Titjdi 



yielding 



L{t) = D,L-^{t), 



L{t) = [iT(0) + 7i?i(t-to)] 



(51) 



(52) 



To determine the decay rate requires a relation between 
fi{t) and L{t) which is provided by invariance of J2: 



J2 = k^{t)L'^+^^{t)C 



(53) 
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with C = f d'^pZfj'{p)T'i{p). Thus, finally, 



(54) 



for [t — to) L''{0)/Di. The generic decay of the mag- 
netic amplitude is predicted by this argument to be deter- 
mined by the scaling exponent which decreases with 
increasing ^ over the range < ^ < Thus, the de- 
cay rate is faster for rougher velocities and slower for 
smoother velocities. It is noteworthy that the decay law 
([54]) is completely independent of the resistivity. 

The above argument does not apply if J2 = 0. In that 
case, one can expect that invariants 



J a,p (^) 



from higher-order zero modes and slow modes of Al* 
(again in the non-gradient sector) will determine the de- 
cay rate . Note, for example, that [d/ dt)J2,i{t) = J2{t), 
so that J2,i becomes invariant if J2 = 0. In ref.[l3| it was 
shown that there are two universality classes in the tur- 
bulent decay of the passive scalar for generic initial data 
with rapidly decaying correlations in space, depending 
upon whether the "Corrsin invariant" Jq from the con- 
stant zero mode is vanishing or nonvanishing. If Jq = 0, 
then there exists a higher-order invariant Ji 7^ 0, asso- 
ciated to the first slow mode in the tower above the 
constant zero mode, which determines the decay. Chaves 
et al. [Tsl showed how this picture emerges from the slow- 
mode expansion of Bernard et al. [l6j and extends to the 
higher-order correlations of the scalar. In the following 
section we shall present a similar treatment of the turbu- 
lent decay of the magnetic field, based on our generalized 
slow-mode expansion in the Appendix. 

There are essential differences, however, between the 
turbulent decay of a passive scalar and of a passive mag- 
netic field. Whereas the scalar field has a finite limit 
as diffusivity k 0, this is not true for the magnetic 
field which, even in the non-dynamo regime, tends to 
accumulate at the resistive scale [2^. As we shall see 
below, the scaling function F*-' (p) grows with p decreas- 
ing through the convective range. Thus, one cannot set 
p = to interpret li^{t) as the magnetic energy. A more 
correct interpretation of the magnetic amplitude is that 
li'^it)/ L'^{t) ~ (|A(t)p), where A is a vector potential 
(2-form) such that B — 8 A. The decay rate of the mag- 
netic energy (|B(i)p) cannot be obtained from purely 
ideal considerations, but requires an explicit matching of 
convective-range solutions with resistive-scale solutions. 
In the following two sections we treat first the ideal, con- 
vective range problem with k ^ 0. 



B. Self-Similar Decay for Initial Data with 
Short-Range Correlations 

Consider any initial 2-point correlation function 
C*-^ (r, 0) of the magnetic field which decreases rapidly 



for large r. We shall demonstrate that the correlation 
C*^(r, t) at much later times exhibits self-similar decay 
and determine the decay law. We use the propagator 
relation Q and the symmetry properties of F to write: 

C={v,t) - j d''pC'''ip,0)FZXp,t\r,0) 

= X' j d''pC^'{p,Q)rU\p,l\v,Q) (55) 

with A = l/{Dit)'^/^ and r ^ r/{Ditf/^. In the last line 
we used the scaling property ((22|) for F. Since C'^^(p, 0) 
decays rapidly for p ^ i(0), we may employ the slow- 
mode expansion dlSl) for {Ditfl^ > L(0). Because 
C'^'^(p, 0) is solenoidal, only the non-gradient zero-modes 
of A^* give a non- vanishing contribution. 
The leading-order term, in general, is 



C^(v,t) ^ {D,t)-(''+~^^^h [ I d''pC^\p^Q)zf,{p) 



rfVC'='(p,0)zr/(p))TyA^ (r,i) 



f(2) 



1/7' 



(56) 



(1) 



In the last line we have used the self-similarity property 
W'(i)(Ar, A'^t) = A-('^+^2'iy(i)(r,t). 

We have also used the fact that W(i) is a real-valued 
function. This will be demonstrated in the following sec- 
tion, where we shall derive the explicit functional form of 
all the self-similar decay solutions. We conclude that, as 
long as J2(0) ^ 0, then the generic magnetic correlation 
C'^(r,t) with short-range initial data is proportional at 
long times to the self-similar decay solution VF^'j*^ (r, i). 

It is important to demonstrate that the above scenario 
is statistically realizable [l^. We shall construct now a 
positive-definite covariance function for which J2(0) ^ 0. 
This will also demonstrate the positive-definiteness of the 
scaling solution , since the dynamics is realizability- 
preserving and the above argument shows that 

lim A''+^2C(Ar, A'^t) = J2(0) • mi)(r,t). 
As a simple example we take, with TV = (o'/V27r)'', 
C''^ (p,0) = Afjd'^k {k^S'^ - k'k^) exp (-^) e'P'^ 

A bit of calculation shows for this example that 
1,(0) = (d-l)(2a2)('^+C2-2)/2^^_^r^'^ + ^2 



X - (C2 + 1)^!^^ , 
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where Sd-i = 27r''/^/r (|) is the hypersurface area of the 

(2) (2) 

unit sphere in d-dimensions and Aj^ , Aj^ are the coef- 
ficients given in the previous section. At generic values 
of d and ^, J2(0) ^ 0. It is noteworthy that J2(0) = 
in the example above precisely at the point of degener- 
acy of zero-modes where (^i = C2 ■ discussed in section 
III C[ this occurs for d > dc = 4.659 at the single value 
^ = < In fact, ^2(0) = at this point for all 

initial data, because there is then a single zero-mode of 

(2) — (2) 

gradient-type satisfying Aj^ = {(2 + l)^Ar ; see 

Of course, whenever J 2(0) — then higher-order terms 
in the slow-mode expansion become dominant and a dif- 
ferent self-similar solution M^(a,p) (r, t) becomes the long- 
time attractor. We shall defer to future work the study 
of this non-generic situation. 

In the remainder of this section we shall make some 
important comments about the generic case J2(0) 0. 
Our first observation is about the property of "quasi- 
equilibrium" . It was shown in [3, that the short- 
distance scaling of the scalar structure function in the 
decay of the passive scalar is identical to the scaling of the 
scalar structure function in a forced steady-state. This 
is the property of turbulence decay traditionally termed 
"quasi-equilibrium." We show here a similar property 
for the turbulent decay of magnetic field, using the slow- 
mode expansion, as in [l^ for the scalar. We use the 
propagator F, its scaling property ((22)) . and the change 
of variables p = \p with A = l/{Dity/'^ to write 

C^{r,t) = J d''pC''{p,0)Fjf,{r,t\p,0) 

= J d''-pC'^'[p/\Q)Fl\{\v^\p,Q). 

We now employ the slow-mode expansion (|42| of F for 
r < {DitY/^ to conclude that 

Ci 



d'p c'^' (pii)'/^P,o) w^if (p,l) 



(57) 



r^(r, t) 



The scaling exponent (^1 and zero-mode are the same 
as found in [26| to determine the short-distance scaling 
of the magnetic correlation function in the forced steady- 
state, which is just the "quasi-equilibrium" property [42l |. 
Since Ci < for all < ^ < we see that C'^{r,t) in- 
creases without bound as r decreases, in agreement with 
our earlier physical discussion. We shall confirm this re- 
sult by an independent argument in the next section. 

A second observation is that the above discussion — 
the demonstration of self-similar decay and quasi- 
equilibrium — carry over directly to the general n-point 
correlation function of the magnetic field. Note that 

t) = (xi, (x2, i) . . . B'- (x„, t)) 



/ 



d''pC^,^^^--^-ip,0)F:X-,ZjAPMr,t) 



where r = (xi,X2, . . .x„) and F„ is the n-body propaga- 
tor. All the symmetries used in the previous argument 
hold for general n, e.g., time-reversal and F„ (p, t|r, 0) = 
Fn (r,0|p, i). Note that due to space homogeneity only 
the separation of the variables matter (the absolute posi- 
tion of each particle is irrelevant) and we can work in the 
(n — l)c?-dimensional separation-of-variables sector. Fn 
then has the scaling property 



Fn{p,0\r,t). 



F„(Ap,0|Ar,A^t) = A-'^f"-!' 

Finally, slow-mode expansions like (|42)) and (|43)) are valid 
for Fn and Fn for all integers n. See jl^ and the Ap- 
pendix for details. The whole analysis thus goes through 
as for n = 2 above and as in [la for the scalar case. 



C. General Self-Similar Decay 

To complement the previous discussion employing the 
slow-mode expansion we shall here determine all possible 
self-similar decay solutions for the magnetic correlation 
function C, following the analysis in 'l4'| for the passive 
scalar. It is convenient to employ the longitudinal corre- 
lation Cl which satisfies the equation ([50]) . We introduce 
the self-similar ansatz 



CLir,t) ^ lr^it)T 



(58) 



Substituting the ansatz ([55)1 into equation ([50)1 for Cl 
arrive at, with p — r/L, 



DiL-^{t) fi{t) 



DiL-^(t) L{ty 



= {d~ i)p«r"(p) + (2e + i)p«-ir'(p) 
H{d-i){d + ^)p^-'np) 

This implies that 

2h{t) 



(59) 



m 

m 
m 



= -aDiL-^{t) 



PD,L-^t). 



(60) 



(61) 



with constants a and (3. We have freedom in choosing 
the value of /3 to fix the length-scale; here we adopt 
13 — 1. The equation for L{t) then becomes identical to 
([5T|) with solution ([52]). Combining ((eO]) and ([61]) yields 
2h{t)lfi{t) = ~aL{t)/L{t) with solution 



^'(t) = [L{t)] 



(62) 



Employing ([SO]) and (pT|) . the equation ([55]) for the 
scaling function T becomes 



[pr'ip)+aTip)] 

= id-i)p^r"ip) + i2^ + d^-i)pT'{p) 
+^{d~i){d + or{p) 



(63) 
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Making the substitution x = —p~'/j{d— 1) yields 

2e 



7^ + 7 + 



d-1 



7^x 



xY^x 



1 



0. 



An equation of this form can be solved by the Frobenius 
method (e.g. see [l^], Sec. 4. 2). According to the general 
theory, there are two independent solutions of the form 
r(a::) — x^^{x), where 6 is a root of the indicial equation 



7 



1) + M + 7 + 



d-i 



If the two roots are distinct and do not differ by an in- 
teger, then the two functions are both analytic, given 
by convergent power series. Otherwise, only one solu- 
tion must be analytic and the second may be an analytic 
function plus Clnx times the first. In our case, it is easy 
to check that the roots of the indicial equation are just 
given by & = Ci/7j '^2/7 in terms of the scaling exponents 
of the zero-modes of A^. If we substitute T = x'-^/^'f into 
the equation for F, we obtain the Kummer equation (28j 



x<S>xx + (c - x)^x - a$ = 0. (64) 



with 



7 7 \ d — 1 

Both independent solutions can be obtained from this 
equation. The first is the Kummer function $(a, c; x), an 
entire function given by the power series 



<I>(a, c; x) 



{a)n x"- 



(65) 



with {a)n = a(a -f 1) . . . (a + n — 1). The other is the 
Kummer function of the second kind, ^'(a, c; x), which is 
defined by a suitable linear combination of ^(a, c; x) and 



l,2-c;x). See [28|, 6.5.6. It is not hard to 



check that this second term corresponds to the root h — 
i-oiH_pi the indicial equation. However, we can argue as 
in [201 that matching the solutions in the convective range 
with those in dissipation range permits only the regular 
zero mode as an admissible physical solution. Thus, we 
obtain Y{x) = x'^i/''$(a, c; x). 

This result can be simplified somewhat by appealing 
to the relation 



C2 = Ci + 7 + 



2i 



1' 



which follows by combining Ci +^2 — —d— 2^/ {d— 1) and 
'^2 + C2 = — rf+7. Note that the above relation generalizes 
the result C2 = Ci + 2 for d = 3 found in [13]. With this 
relation we obtain c = (Ci + C2 + d)/'^, so that 



F(p) = p^i$ 



a + Ci C1+C2 



P' 



(d-l)7 



(66) 



All the self-similar solutions of Eg. ([30)) are given by 
the ansatz ([55[) with a scaling function of the form in 
([66[) above and with L{t) and ti{t) given by Eqs. ([52[) 
and (|62p . respectively. Since ^•(0) = 1 is finite, all of 
these self-similar solutions satisfy the condition of "quasi- 
equilibrium," showing the same scaling r^^ for r <C L{t) 
as found in [26| for the forced steady-state. 

There are two distinct types of self-similar decay so- 
lutions corresponding to different choices of the free pa- 
rameter a. When a = C2 + d+pj, for p = 0, 1, 2, . . ., then 
a — c + p with p = 0, 1, 2, . . .. In this case 



<I>(c + p, c; —x) 



^\x)e 



where Lp^^(a;) is the generalized Laguerre polynomial of 
degree p. (See [111, Ch.6). This series of solutions has 
stretched-exponential decay in space. If, for example, we 
take a = C,2'^ d corresponding to p = 0, then we get 



F(p) = exp 



1 



1 7 



The corresponding self-similar decay solution satisfies 
C{\v,X<t) = A-(''+«2)c(r,i). The a = C2+c^ solution thus 
coincides with the self-similar solution W(i)(r, i) which 
was shown in the previous section to describe the long- 
time decay of generic initial data with short-range corre- 
lations. More generally, the solutions with a = C,2~^d+^p 
coincide with the self-similar solutions W(i_p)(r,f) for 
p = 0, 1, 2, . . . which appear in the slow-mode expansion 
p3[) of the adjoint propagator F. 

For any other choice of a 7^ C2 + + 7P with p = 
0, 1, 2, . . . one obtains instead a class of self-similar de- 
cay solutions with power-law decay of correlations at 
large distances. This follows from the asymptotic re- 
lation <^ia,c\~x) ^ tJ-^^cc"" for Re a; +00, if 

a ^ c -|- p, p = 0, 1, 2, . . .(|23|, 6.13.1). Using the above 
relation together with ([M[l . ([55[) . and ([S^ gives for any 
self-similar solution with a 7^ C2 + d + ^p, p = 0,1,2, ... , 



r > L{t), 



where A is a time-independent constant. This result is 
usually called the "permanence of the large-scale eddies" 
in the turbulence literature. Note that for initial data 
with such power-law decay of correlations, the relation 
between fi{t) and L{t) that determines the decay rate is 
obtained from this permanence, as fi^{t) ~ A[L{t)]~°' , in 
agreement with (l62t . See [l3| for more discussion. 



D. Decay Law of the Magnetic Energy 

We are now ready to discuss the decay law for the 
magnetic energy: 

£;(t) = i(|B(i)P)=trC(0,i). 
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Under the assumption of isotropic statistics made here, 
E{t) = (d/2)CL(0, t). Clearly, in order to evaluate this 
expression at r = 0, we must consider the matching of 
our convective range solution to the resistive scales. We 
may do this heuristically, as follows. We assume that, to 
leading order, 

E{t) {d/2)CL{£.,t). 

We then estimate the correlation function on the right 
by matching with the convective-range expression 



CL{r,t) ^ Coh\t) 



coim] 



for r ^ L{t) and some positive constant Cq. This yields 

for Ci = {d/2)Co. Although the energy magnitude in- 
creases as resistivity is lowered, the decay rate is in- 
dependent of resistivity. Since Ci < 0; we see that 
the decay of magnetic energy E{t) (x {t — iQ)-("+Ci)/T 
is always slower than the decay of magnetic amplitude 
n'^it) (X {t — ta)~°'^^ . For example, for the generic case 
with a = C2 + wc obtain that 

E{t) cx {t-to)-", 

with c= (Ci +C2 +c^)/7- 

The above argument is basically correct, but not fully 
rigorous. It seems worthwhile to give a more systematic 
derivation, since the time-dependence of magnetic energy 
is crucial to the issue of whether dynamo action is present 
or not. We employ a standard method of matched asymp- 
totic expansions (see, for instance. Chap. V of Ref. [29j). 
The equation obeyed by Cl for k > is 



dtCL = [{d-iy^dlCL + i^^ + d" -I) 



drCL 



-2k 



1 



d^CL + {d+l)-drCL 

r 



.(67) 



See 0. Substituting the self-similar ansatz ([55]) we get 



K + pTp] + [{d - l)p«rpp + (2^ + d^- ijp^-'Tp 
+e(d-l)(d + Op«"'r] +e« 



' (68) 

with p = f/L{t) and e = £^/L{t). In the outer range 
where p — 0(1), the dominant balance in the equation 
(j68p is between the first term from the time-derivative, 
which acts like a forcing, and the second term from the 
turbulent advection. The third term may be neglected for 
small e , yielding the leading-order equation for the outer 
solution. This is the same equation which was examined 
in the preceding section UlI CI where all self-similar solu- 
tions were found. Thus the outer solutions ro„t(p) are 
given by the formula (j66p for any choice of a and multi- 
plied by an arbitrary constant Cout ■ These solutions must 



0. 



now be matched to an appropriate inner solution in the 
resistive range. 

We introduce the inner variable a — r/^K = p/ e in Eq. 
to obtain 



[ar -f aT„] + [{d - l)a^T„„ + (2^ + d^ - l)ai-'^T„ 

Hid ~ i){d + O'^^^T] + [r,, + {d + i)ir,] = 0. 

(69) 

The dominant balance in (j69p is between the second term 
from the turbulent advection and the third term from the 
molecular resistivity. To leading order we can disregard 
the first term proportional to to get 



a^r,, + {d+ 1)C7T, + a« [{d - l)a^T,„ 

+ (2^ + d^- l)aT, + i{d - l){d + or] = 0. 



(70) 



Making the change of variables y = —{d— \)a^ reduces 
this to a hypergeometric equation ([11], Ch.II): 

y(l - y)Tyy + [c, - (a* -I- 6, + l)y\ Ty - a^bj" = (71) 

where 

a, + 6, = - I Y + d 1 , a^b^ = = ^ . (72) 

Up to an overall multiplicative constant, there is a unique 
solution F{a^:,b^]c^,;y) of the above equation which is 
analytic in the region arg(l — y) < tt of the complex y- 
plane. This hypergeometric function is given for \y\ < 1 
by the absolutely convergent power series. 



F{a^,b^;c^;y) = ^ 



71=0 



{a*)nib*)n 

(c*)„ n! ' 



(73) 



if 7^ 0,-1,-2,.... The other independent solution, 
y^-^'F{a, + 1 - c,.,b, + 1 ~ c,;2 - c,;y) [HH, 2.9(17)], 
is singular at y = and must be discarded. Because of 
the symmetry F{a^,b^;c^\y) = F{b^,a^\c^;y) we have 
freedom in choosing a* and 6* . Combining the equations 
in (|72p yields a quadratic equation for a.^, 

fi* ^ T — 7 + " a* H T — = 

i\d-l J ^ 



and an identical equation for b.^, . It is easy to check that 
the roots are just — Ci/C ^^^id —i^iji,, where Ci7^2 are 
the scaling exponents foimd in section Hi CI We choose 
a* — — Ci/'f and 6* = —ujij^^. Thus, we obtain 



r„(a) = a.„F(-|,-^^;i±^;-(d-l)a« 



for the inner solution, with an arbitrary constant Cin- 
This solution gives the complete description in the re- 
sistive range, e.g. implying a magnetic energy spectrum 
E(k) cx for IJi > 1. 

To match this solution to the outer solution, we must 
find its asymptotic behavior for cr ^ 1. This is given by 
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F{a^,K;c^]y) - r(bl)r(cl-al) (-?/)~°'' as Rey ^ -oo, 
for a, < 6*, a, ^ c, + p with p = 0, 1, 2, . . . [see psf . 
2.1.4(17)] to be 



r(c,)r(&. 



r(6*)r(c* — a*) 



(d-i) 



for (7 3> 1. This is the same power-law as Toutip) ^ 
Coutp'^^ for p ^ 1. Equating the inner and outer solutions 
Tin (c) = Tout {p) in the overlap region e ^ p <C 1 yields 
the relationship 



Cin — 



r(6,)r(c* 



r(c*)r(6* 



Notice that the first factor is a numerical constant B{£) 
satisfying B{Q) = d — 1 and -B(^*) = 0, and varying 
smoothly between those limits. 

Finally, we obtain the magnetic energy from E{t) = 
(d/2)^2(i)r,„(0) = {d/2)fi'^(t)Cin which, with e 
and = [L(t)]-", gives 

E{t) = Ci[L(t)]-'"+'^i)£^i 

for Ci = {d/2)B{£) ■ Gout- This differs from the previous 
heuristic estimate only by a constant factor. 



E. Magnetic Induction and Dynamo Order 
Parameter 



The above arguments are reminiscent of our discussion 
in subsection III C[ where we emphasized the importance 
of considering the correlations between line-vectors ad- 
vected to the same point, in order to distinguish between 
dynamo and non-dynamo regimes. In fact, the two sub- 
jects are intimately related. As we now show, the "or- 
der parameter" 'R-{t) that we considered in ([SD]) can be 
interpreted as the energy of a certain self-similar decay 
solution C(o) corresponding to an initial condition which 
is a random, statistically isotropic but spatially uniform 
magnetic field. Such a random magnetic field has a co- 
variance of the form 

for a positive real number A. A constant correlation such 
as above would be invariant for an advected scalar, but it 
is not for a magnetic field. There is a well-known phys- 
ical phenomenon of "shredding" ^] or "induction" |12l | 
of a constant magnetic field due to the stretching term 
(B- V)u in the evolution equation. Thus, an initially con- 
stant magnetic field will develop very fine-scale structure 
by turbulent induction and may — in principle — act as a 
seed field for kinematic magnetic dynamo. 

The correlation at later times with the above initial 
condition is provided by ([H]), which yields 



For the limit k — > in the KK model, the scaling relation 
([22ll for F then implies that 



(r,t). 



Thus, C(o) is a self-similar solution of 9tC(o) = A^C(o). It 
is clearly the self-similar solution with parameter a = 
in our general classification. On the other hand, if we 
take 1/d then also 



= \jd''pFg{v,t\p,Q) 

= {5t{t)5e^{t))^. 



The latter expression denotes the correlation of line- 
vectors which started as the same random unit vector 
at time 0, at any pair of points, which ended up at time t 
at points displaced by r. We should emphasize that this 
result is valid for any divergence-free advecting velocity 
field and thus applies as well to incompressible fluid tur- 
bulence. It immediately follows by summing over i = j 
and setting r — that 

7^(^) = 2i;(o)W 

where i?(o) {t) is the energy of the solution C(o) . 

Our analysis in the previous section can be applied 
to describe the behavior of C(o)- The formula CL{r,t) — 
T{r/L{t)) holds using the analytic expression ([66]) for F 
with a = and L{t) = [Ditf/^, vahd for all r » 
It is more interesting to consider various asymptotic be- 
haviors. The "permanence of large eddies" implies that 

CL{r,t)^A, r:$>L{t). 

In the convective range 



CLir,t)^A 



L{t). 



Finally, for r — > and long times, 

which grows with decreasing k or increasing t, but only 
as a modest power law. 

This result may be interpreted in terms of material-line 
correlations by setting A — 1/d : 

n{t) - {5i{t)-6t'{t))a cx li'{D^t)\^'\'\ (74) 

which implies that this quantity grows slowly with time. 
It would be of great interest to determine the time- 
dependence also in the dynamo regime. If the leading 
eigenfunction of A^* satisfies / d'^r£ki{v) = 0, then 'R{t) 
need not grow exponentially. Note, for example, that 
the space-integral of the dual eigenfunction £*-' (r) does 
vanish, so the issue is not straightforward. 
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IV. FINAL DISCUSSION 

Our work leads to several important conclusions re- 
garding the small-scale turbulent kinematic dynamo. 



A. Breakdown of Flux Freezing and Dynamo 

In order to understand the turbulent dynamo process a 
crucial fact is that magnetic field lines are not frozen into 
the plasma flow, even in the zero-resistance limit k — > 0. 
Flux-freezing would imply that only a single field line is 
advected into each space point from the field configura- 
tion at an earlier time. In fact, infinitely-many field lines 
are carried to each point by a combination of fluid ad- 
vection and resistive diffusion [l^, [l3] ■ In the Kraichnan 
velocity ensemble, the probability for two line elements 
to arrive at the same point at time t starting from points 
separated by r at time is P{0, t\r, 0) cx exp{—r'^/^Dit) 
in the limit k — > and does not degenerate into a delta 
function S'^{r) This is a manifestation of the phe- 

nomenon of "spontaneous stochasticity" first pointed out 
by Bernard et al. [l^ , which is due to the explosive sep- 
aration of pairs of fluid particles undergoing turbulent 
Richardson diffusion. It was argued in ^iu\ that this be- 
havior as K — > holds in general for a turbulent plasma 
with a rough velocity fleld, so that Alfven's theorem on 
flux-conservation remains as a stochastic law only. 

The breakdown of flux-freezing in the case of rough 
velocity fields renders the turbulent kinematic dynamo 
an even more subtle problem than the laminar (or large 
Prandtl number) kinematic dynamo (for the latter, see 
e.g. refs. [3li|3^,l33|,|34,|35|-) For the very smooth veloc- 
ities considered there (^ = 2), Alfven's theorem holds in 
its usual form in the limit k 0. However, for rougher 
velocities with rugosity exponent anywhere in the range 
< ^ < 2, an infinite number of field lines enters each 
point even in the zero-resistance limit. The resultant 
magnetic field is the resistive average over the field vec- 
tors of all of the individual lines. We have shown in this 
work that the presence of small-scale kinematic dynamo 
effect depends upon the existence of sufficient angular 
correlation between the individual field vectors. Thus, 
dynamo action occurs in the KK model for smoother ve- 
locities with ^* < ^ < 2 but not for rougher velocities 
with < ^ < This is true despite the fact that the 
stretching rate of individual fields lines is much greater 
for ^ smaller. 

In section HIC] we defined 7^fc£(r,^) = i^*^(0, t|r, 0), 
which measures the correlation between line-elements 
6£k(t) and S£i{t) at the same point at time t which 
started out as unit vectors e^, eg at distinct points sep- 
arated by r at time 0. We found there that, in the non- 
dynamo regime of the KK model with < ^ < C*, this 
quantity scales as 



for iiy,£K r <^ (Dit)^/'^. Here Z is an appropriate 
zero-mode of M* scaling as Zke{r) oc r'', with ~d < C < 
0. [Note that the factor arises from VF*'(0,1) in the 
slow-mode expansion.] This correlation decays only as 
a power for r increasing through the inertial-convective 
range, implying that line vectors initially separated by 
distances ~ (£)ii)^/^ contribute substantially to the final 
average. The correlation does not vanish as k — > but, in 
fact, increases as a moderate power of i^, demonstrating 
that infinitely-many field lines continue to contribute in 
that limit. The result is, however, a correlation TZm{v, t) 
slowly decaying in time. On the other hand, the lengths 
of the individual line-elements {S'^£k{t)), {5'^£i{t)) grow 
exponentially in time as in l|37p with rate A oc v/£l, = 
\/t^. The result is that 



0, 



(75) 



7^fe,(r,t) ~ Cii}{Dit)-^ZM{v), 



exponentially rapidly either as t ^ oo or as k ^ with 
u < PtcK. We conclude that the dynamo fails for a very 
rough velocity field because advected line-vectors arrive 
at the same point with insufficient angular correlation. 
Although individual field-lines are stretched to an in- 
credible degree, resistive averaging of nearly uncorrelated 
lines leads to almost complete cancellation. 

The situation is qualitatively different in the dynamo 
range with smoother velocities (5* < ^ < 2). In that 
case, we have from (|49|) that 

7^fc4r,t)cxe^°*ffc£(r), 

where Eq is the dynamo growth rate. Since Eq cx l/t^ <C 
A (X 1/ty, for A in ([55]) . it is still true that the angular 
correlations (j75p decay exponentially either as t — > oo or 
as K — > with Pr small enough. However, the decay 
exponent is reduced by a finite amount. Enough cor- 
relations remain between line-elements entering a point 
that the net magnetic field after resistive averaging can 
profit from stretching of individual lines and exponential 
growth of magnetic energy ensues. 



B. Hydrodynamic and MHD turbulence 

Much of the formalism of this paper carries over to 
the problem of kinematic dynamo for a weak seed mag- 
netic field in hydrodynamic turbulence. The propagators 

F^^i{p,t\r,0) — Fll{r,0\p,t) give a fundamental descrip- 
tion of the kinematic dynamo for any incompressible ad- 
vecting fiow. All of the results of section III Al apply in 
general, in particular equations Q, pT|) and p^ . and 
also the relationship in section llll El between magnetic in- 
duction and line- vector correlations. Any further simpli- 
fications from space-homogeneity and time-stationarity 
also apply where appropriate. On the other hand, some 
features of the KK model are quite special and do not 
apply more generally. The self-similarity property 
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of the propagators F and F does not carry over to hy- 
drodynamic turbulence, because of small-scale intermit- 
tency of the advecting velocity field. Also, the statistics 
of forward and backward Lagrangian trajectories are not 
identical in hydrodynamic turbulence [36[. Thus, rela- 
tions such as (fTHl) which depend upon time-reversal sym- 
metry do not apply to the real turbulent dynamo. Lastly, 
the time-evolution of the propagators F and F is in gen- 
eral non-Markovian and thus the simple diffusion equa- 
tions such as (P(H) . (PT|) do not apply. One of us (G.E.) 
is currently carrying out a numerical evaluation of these 
propagators for hydrodynamic turbulence, which will be 
reported elsewhere. 

We expect that many of the ideas of this work will ap- 
ply even to nonlinear MHD turbulence and dynamo ef- 
fect there. A stochastic form of flux- freezing and Alfven's 
theorem holds also for non-ideal (viscous and resistive) 
MHD [l3l. We expect these conservation laws to re- 
main stochastic in the limit k — > 0, — > with Pr 
fixed js^]. However, there will be nontrivial differences 
from the kinematic problem studied here, due to backre- 
action of the magnetic field on the plasma flow via the 
Lorentz force. For example, in comparison with hydro- 
dynamic turbulence, 2-particle relative diffusion in MHD 
turbulence is observed to be suppressed in the direction 
transverse to the local magnetic field [ST*]. In principle, 
however, one can account for all such nonlinear effects 
by the presence of a second stochastic conservation law, 
a modified Kelvin theorem [l3, [H, [3p| . We believe that 
"spontaneous stochasticity" and the implied stochasticity 
of magnetic-line motion and flux-conservation must play 
a central role in the understanding of MHD turbulence, 
dynamo and reconnection. 
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V. APPENDIX: SLOW MODE EXPANSION 
FOR NON-HERMITIAN EVOLUTION 

Unlike for the passive scalar, the n-body evolution op- 
erators A4„ for the passive magnetic field in the Kraich- 
nan model, are no longer even formally Hermitian, with 
Al* ^ Mn- Nevertheless, certain important properties 
of the scalar evolution operators remain true for A1„ 
and Mn- these are homogeneous of degree —7, reality- 
preserving, and — in the non-dynamo regime — having ab- 
solutely continuous spectrum on the negative real axis. 
As we shall show in the following, the above properties 
together with assumed analyticity conditions allow us 
to generalize the zero-mode and slow-mode expansions 
derived in [iGj for the Hermitian case to pairs of non- 
Hermitian operators Ai and M*. Although our intended 
application is to the Kazantsev-Kraichnan kinematic dy- 
namo model, we shall carry out the derivation in an ab- 
stract, general setting. We shall employ the properties of 
Ai and A4* given above and, also, other properties that 
will be stated explicitly below. The entire argument is 
modelled very closely after that in [l^, with just a few 
important differences that are stressed below. 

A. The Zero- Mode Expansion 

We shall assume that the operators M and Ai* act on 
i^(M''). The dimension d need not be identified with the 
physical space dimension, as in the main text of the pa- 
per. (E.g. if ds is the space dimension, then the n-body 
operators in the Kraichnan model act on L'^{M.'^) with 
d — nds or with d = (n— l)d in the translation-invariant 
sector.) Define Green's functions for the operators by 

-A^xG(x,y) = -Al*G(x,y) =(5'^(x-y), 
-Al*G(x,y) = -XyG(x,y)=^'^(x-y), ^'"^ 

where the subscript (x or y) indicates on which variable 
the operator acts. Note that these Green's functions are 
both real-valued and, of course, G(x, y) — G(y,x). 

Our aim is to derive the following short-distance 
asymptotic expansion for G : 

G(x/L,y)^^L-';»/,(x)[g,(y)]*, L » 1, (77) 

a 

where * here denotes complex-conjugation. The function 
fa is a regular zero-mode of Ai with scaling dimension 
Ca, while is a singular zero-mode of Ai* with scaling 
dimension uja — —d + 7 — Ca- What dominates in the 
expansion (|77p is the contributing zero-mode whose scal- 
ing exponent Ca has the smallest real part. Thus, we 
label the exponents according to the magnitude of their 
real part, so that Re (a > ReCa' and Rewa < RecJa' for 
a > a' . We derive also a similar expansion for the adjoint 
Green's function 

G(x/L,y)^^L-^"7jx)[5,(y)]*, L » 1, (78) 

a 
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where now is a regular zero-mode of A^* with scahng 
dimension while ga is a singular zero-mode of M with 
scahng dimension uja = —d + j — (^. We thus see that 
the homogeneous zero- modes of the operators A4 and A^* 
come in pairs, {faT9a) and {fa,'ga)^ with related scaling 
exponents. 

Following [rG*! , we employ the Mellin transform, which 
is a unitary transformation between the spaces L^(M'') 
and L^{Rea = -d/2) ® L^^^-d-i) ^^^^^ 



(79) 



at complex values Ca, a = 1, 2, ... in the strip. By mov- 
ing the integration contour in (j84p further and further 
to the right, one picks up successive pole contributions. 
This implies that Green's function for large L satisfies: 



/(x)K^/(a,A)= / A— V(Ai)dA 



with the inverse transform, for R = Ixl 



with the function 

Z„(x,y) = [i?(x)]^"Z,(x,y)[i?(y)]- 



■<i-l-7-Ca 



R''f{a,±)da. 



(80) 



which is homogeneous of degree Ca in x and of degree 
—d -|- 7 — Ca in y. From the definition of the Green's 
function, using A^x = i~^A^x' with x' = x/i. 



-L-'^S'^ (f 



The inner product on L^{Rea = -rf/2) (g) i^(5''"^) is 



-^xG($,y) 
E,i^-'^"Xx^a(x,y) 



f^9) = 



dti;(x) 



da 



Re a—- 



f{a,±) .g(cr,x) 
(81) 



from which we get A^x^a(x, y) = for points off the 
diagonal. Likewise, 



However, it is more convenient to write this as 



f,9 



da 



Re a- 



/(-cr*-d,x)J 5(cr,x) 
(82) 



Although 



/(^,x) 



on the line 



/(-a*-d,x) 

Re (7 = — rf/2, the second expression is analytic in a when 
/((T, x) is analytic. This form of the inner product allows 
one to shift integration contours in the complex cr-plane. 
A key role in the analysis is played by the operator 

R'^/'^MR''/^ 



(83) 



which is homogeneous of degree zero. Since it thus com- 
mutes with the self-adjoint generator D = j (x- Vx -\- f ) 
of dilatations, it is partially diagonalized under the Mellin 
transform: 

(AA/)~(a,x)=Ar(fT)/>,x), 

where N{a) for each a is an operator on Li^ {S'''~'^) . Using 
M.~^ = i?'^/^A/'~^i?^/^, one straightforwardly derives the 
following fundamental identity for the Green's function 
G(x,y) --A^-i(x,y): 



from which we get A^*Za(x,y) = for points off the 
diagonal. We finally conclude that Za(x, y) for fixed y 
is a homogeneous zero-mode of A^x of degree and for 
fixed X is a homogeneous zero-mode of A4* of degree 
—d -l- 7 — Ca- If we assume that zero- modes of a given 
scaling exponent are non-degenerate, as will generically 
be true, then we can write 

Z,(x,y)=/,(x)[gjy)]*, 

where fa is the unique scaling zero-mode of M. with ex- 
ponent Ca and is the scaling zero-mode of A^* with 
exponent ZUa = — d -I- 7 — Ca- We have used here the fact 
that Al* is reality-preserving. This yields ([77|) . The ex- 
pansion (|78p for G is derived by an identical argument. 

Although we shall not employ the corresponding large- 
distance expansion in this work, we make here a few re- 
marks about it. Under the Mellin transform the adjoint 
of has the kernel 



G(x,y) 



Rea—- 



7 . . 

2;x,y 



(,c^;x,yj 



AA-i(-a*-d;y,x) 



X [i?(y)]-'^+^-° . (84) 

See [lB|. We note that the shifts in cr arise because R^/'^ 
acts as a translation by —7/2 under the Mellin transform. 
The above identity is the key to deriving the zero-mode 
expansion for G. 

The main hypothesis is that the operator function 
M^^{<7) is meromorphic in a wide vertical strip around 
the line Rccr = — rf/2, whose only singularities are poles 



This last relation reveals the important fact that if 

N^^{a) has a pole at Ca then M* (cr) has a pole at 
TDa = —d + 7 — Ca- Indeed we have: 



(a - ^;x,y) 



AA-i(-d+i-a*;y,x) 



-N- 



7 



:,x,y 



^a(x,y) 

- Ca 



(-<i+7-<T*)-C< 

-1 TT 



-Uy)[gMr\ 
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At this pole with KeZJa < —d/2 the role of the regu- 
lar and singular zero-modes in the residue is reversed. 
By pushing the integration contour in (|84p further and 
further to the left one can thus derive a large-distance 
expansion for G. This can also be directly obtained from 
the short-distance expansion for G, as follows: 

G(Lx,y) - G(ix,Lf) 
= L^-'G{^,i) 

= L^-'^^L-C[/4y)]*5jx) 

a 

a 

for L ^ 1. Of course, a similar expansion holds for G. 



The Slow Mode Expansion: Elementary 
Arguments 



Define the heat-kernels 

P(x, i|xo,io) = (x|e 
P(x, i|xo,io) = (x|e 



{t-tQ)M 



xo), 



(85) 



so that, obviously, -P(x, i|xo, to) — P(xo, i|x, to). We have 
the relations 



G(x,y) = / dtP(x,<|y,0) 
Jo 



and 



G(x,y) 



dtPi^,t\y,0). 



(86) 



(87) 



Given the validity of the zero-mode expansions for G and 
G one should expect that related expansions hold for P 
and P. We shall show that this is indeed true, with the 
asymptotic expansion analogous to (|77p for L ^ 1 : 



^^(^,^1x0,0 



E 

a,p>0 



i-(^°+^^Va,p(x)[5a,p(xo,t)]' 



(88) 

Here /a,p are the tower of regular slow modes of M, satis- 
fying -Mfa,p = fa,p-i and fa,o = fa- Also, g^^ p are solu- 
tions of 9tg„ p(x, t) = A^*g^ p(x, t) with initial conditions 
5a,-i(x,0) = g„(x) and g„_p+i(x,0) =^ ->l*5„_p(x, 0). 
They satisfy the scaling relations g^p{\x,\''t) = 
X^c^-(p+'^)tg^^(^x,t). Note that the dominant contribu- 
tion in (|88[) will generally come from the tower with min- 
imum Rc(C(j) and from the first (zero-mode) term p — 0. 
There is an analogous expansion for P with L 3> 1: 



P(^,t|xo,0 



^ L-(C.+7rt/^_^(x)[g,.p(xo,t)]*, 

a,p>0 

(89) 



with the roles of the operators Ai and Ai* reversed. 

We shall derive the above expansions in this section 
and the next. Here we proceed by assuming that a gen- 
eral expansion exists for L » 1 of the form 



Pg,t|x',0) -^L-''"/„(x)[5jx',t)]' 



(90) 



We shall then identify the form this expansion must take. 
In the following section we establish from a more funda- 
mental point of view the existence of such an expansion. 
First we substitute ([M)) into 

atP(x, t|x',0) = Ai;,P(x, i|x',0) 
= AlxP(x,t|x',0) , 



obtaining 



J^i-^^/^x) [9t5.(x',i)]* 

= ^L-''"/„(x) [X;,5jx',t)]* 

a 

= J2L-P-+''MM^)[9M,t)r 



(91) 



We see that whenever the asymptotic series contains a 
term proportional to /q(x) with scaling exponent pa it 
must also contain a term A^x./q(x) with exponent Pa — 1, 
and then a term A1J/q(x) with exponent pa — 27, and 
so on. This cannot continue indefinitely, since, other- 
wise, there would be successively more and more diver- 
gent terms for L 1. The only way that this sequence 
can terminate is if, eventually, 

Al£+Va(x) - 

for some integer p. In that case, we see that fa = 
i—Mx)^fa = fa,p for some homogeneous zero-mode fa, 
and the expansion (l90|) contains the whole tower above 
that zero mode. All such towers associated to regular 
zero modes must appear because the condition (|86p to- 
gether with the zero-mode expansion for G implies that 

r POO 

^L-^^/^x) / dtgj^\t) 

-^L-^Va(x)[5.(x')]*. (92) 

a 

The expansion (|5D|) thus must have precisely the form of 
equation (jSSp and we must only establish the properties 
of g^ = ga,p- We note from ([9T|l that 

dtVa^p = M*ga,p = -5a,p+l- 

Also ((92)) implies that (away from the origin x' = 0) 

/•oo 

5a,p-i(x',0) = - / dtdtga,p-i{x.',t) 
Jo 

/>oo 

= / dtg^^j,{K',t) = 
Jo 
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for p = 1,2,3..., whereas gg^_i{x',0) = 5a(x'), the 
singular zero- mode of M*. Finally, the scaling prop- 
erties of p follow from the scaling property of 'g^ 
and of the propagator, i.e., g^(Ax) = X'^''g^{x) and 



(Ax, Ay) = A ''e (x, y), respectively. The ex- 
pansion ([89|) for P is derived by an identical argument. 



The Slow Mode Expansion: Fundamental 
Derivation 



D. It is convenient to introduce instead the self-adjoint 
operators K{a) = U-^{a)Ko{a)U{a). Thus, 



X 



U{(T)[-id„+iK{a)]U-\a), 



(100) 



X* — >U{a)[-fd„-iK{(j)]U-^{a). (101) 

We now introduce the operators L{a) on L'^{S'^^^) sat- 
isfying 



We shall now demonstrate the existence of the expan- 
sion ([50]) and verify by an independent argument its gen- 
eral properties discussed above. A key fact that we use 
is that the operators M and both have spectrum 
absolutely continuous over the negative real axis. This 
assumption explicitly rules out kinematic dynamo effect 
due to point spectrum on the positive real axis. Because 
of this assumed property, we may define 



X = log(-7W), X* = log(-M*) 



(93) 



where the branch of the natural logarithm log(z) is de- 
fined with cut along the negative real axis. Furthermore, 
because M and Ai are homogeneous of degree —7, the 
operators X and X* both satisfy the Heisenberg commu- 
tation relations 



[D,X] = [D,X*]^i^I, 



(94) 



where D is the self-adjoint generator of dilatations. 
We may decompose X, X* into Hermitian and skew- 
Hermitian parts, as 

X = H + iK, X* = H-iK, (95) 

where H, K are both Hermitian. In that case, we see that 

[D, H] = i-fl, [D, K] ^ 0. (96) 

We can now follow the arguments in Ref. 1^ to infer 
that under the unitary Mellin transform 



D 



1 / d 

ir+2 



H 



f/(a)79,[/-i(a). 



K 



Ka{cj) = U{a)K{a)U-\a), 



(97) 



(98) 



(99) 



where Re{a) = — |. As in Ref. [l^, the operators U{a) 
in Eq. ([98|) are unitary operators on LF'{S'^~^) and the 
result in Eq. (j98p follows from the Stone-von Neumann 
theorem on uniqueness of representations of the Heisen- 
berg algebra. The operators ifo(<T) in (|99p are self-adjoint 
operators on L^(S''^~^) and the result ([M|) is a conse- 
quence of the second half of ([M)) — commutativity of D 
and K — so that K leaves invariant the eigenspaces of 



^ L{a)^~iK{a)L{a), 1(0) = I, (102) 
da 



-f—L-\a)=L-\a)iK{a), L-\0)=I. (103) 
da 

The operators L{a) and L~^{a) can be defined explicitly 



by ordered exponentials along the line a — 



L{a) = 




and 



L-\a) 




-l.f^d.'K{-i+^.') 



It follows that 

-fd, + iK{a) ^ L{a)-idaL~\a), 

jd, ^ tK{a) = L*-\a)^d„L*{a). 



if > 0, 

if 1/ < 0. 

(104) 

if > 0, 

if < 0. 
(105) 



(106) 



(107) 



Finally, combining (fT06| . (fT07| with (fTOO]) . (fTOT|) . we 
obtain the mappings under the Mellin transform 



X 



V{a)jd,V-\a), 



X* 



V*-\ahd,V*{a). 



(108) 



(109) 



with 



V{a) ^U{a)L{a), V*{a) = L*{a)U~\a). (110) 

This is the main result that we require. 

The rest of the derivation of the slow mode expansion 
follows the argument of Ref. [l^, assuming that V{a) ex- 
tends to a meromorphic operator-valued function of a. 
We shall sketch here the main points. Note first that we 
can exponentiate the relations (|108p . (|109p to obtain 

-M = VR-''V~^, -M* = V*-^R-''V*, (111) 

where we have defined the operators V and V* by 

(F/)-(a,i) = V{a)f{a,i), {V*fr{a,i) = V*{a)f{a,i), 
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which are mutual adjoints. From the definition ^ 
i?-7/2_/V/t- 1^-7/2 and PTT|) we see that 

which under MeUin transform becomes 



Of course, we have also 



-M*-'(a-l)=V*-\a)V*{a~^), 



(112) 



(113) 



by an identical argument. 

One immediate consequence of (|112p is that poles of 
(cr — ^) can arise only from poles of V{(j) or zeroes 
of V{a — 7). Our main assumption will be that all of the 
poles of V{a) lie in the half-plane Rccr > —d/2 and all 
of its zeroes lie in the half-plane Recr < —d/2. Because 
of the adjoint relation 



V* 



* -1 



(114) 



we see that V* "^(cr) then enjoys the same property, with 
the poles of V{a) corresponding to zeroes of V* ~^ (cr) and 
the zeroes of V{(j) corresponding to poles of V^*^^(ct). 
The assumption on V implies that all of the poles of 

poles of V{a) with the form 
1 



^) for Recr > — (i/2 + 7/2 must arise from 



V{a) 



l/a)(5al^K-7), 



(115) 



in order to reproduce the known poles oi Af ^ (c — 
On the other hand, we can rewrite (|112p as 



V{a + -fp) = -7V"^ cr + 7 U 



y(a + 7(p-l)), 



Finally, we exponentiate one more time relations (jllip 
to obtain 



tM 



e-- = T/e-*«"V-\ e*-^' =l-*-ie-*^"V*. (116) 
The first of these, under Mellin transform, gives 



(e--^)(x/i) = 

1 f da 



7 JRcCT=-d/2 
X 



27ri 



/ dco(y)y(a;x,y) 



Rco-' = -d/2-0 ^Tl"'* V 7 



Pushing the cr-integration contour further and further to 
the right gives the expansion for L ^ 1 : 

(e*^<^)(x/L) - ^L---^''/,,,(x)(g,^p(t),^), 



with a suitable definition oi'g^p{t) See [lg| for more de- 
tails. The above is just an integrated form of the slow- 
mode expansion (I88p for P. The slow-mode expansion 
(155)) for P follows by an identical argument, in which 
the various terms arise from the poles of V* ^^(cr) in the 
half-plane Recr > —d/2. 

There are similar large-distance expansions for P and 
P, in which enter the "tunnels" of singular slow modes. 
The terms in these expansions arise from the zeroes of 
V{a') and V*-^{a') in the half-plane Recr' < -d/2 by 
moving cr'-integration contours to the left in a formula 
similar to the above. The reader may work out details. 
We shall just note here that the pole (jllSp of V^(ct) implies 
via the relation 



V*{a) = [V{-a* -d)]* 



the result 



F*(a-7) = 



-V*{u)a)\g,){fa\ 



(117) 



for p = 1, 2, .... Let us assume that none of the poles ^^"^^ ^'^y" ^* ^ 

of ('^ ~ 2) occur at points in the complex cr-plane 
with real parts differing by integer multiples of 7. This 
will hold generically. In that case, (o-q -I- 7 (p — 5)) 
is a regular operator for all p = 1,2,... and we may use 
the above relation to infer inductively a series of poles 



1 



cr - Ca - 7P 

for each a = 1, 2, . . . with 



\fa,p){9a\V{<Ja-l), 



fa,p-l 



71P-2 



for p = 1, 2, . . . and fafl — fa- It is not difficult to check 
that this coincides with the definition of fa.p given earlier. 



— 7. This zero 

and the "tunnel" of zeroes beneath it give rise to the 
terms in the large-distance expansion of P. 

In our discussion throughout the appendix, we have 
assumed that all the regular zero modes of A4 and M* 
have scaling exponents a with Recr > —d/2 + 7/2 and 
all the singular zero modes have scaling exponents with 
Recr < —d/2 — 7/2. This is true in the KK model only 
for d > 6 and for ^ not too large. For d < 4 the two 
primary "singular modes" have exponents , tUi > —d/2 
for all ^ and, for sufficiently large ^, these exponents even 
cross and become larger than Ci j Ci j respectively! For 
d > 6 it still happens that Ci < —d/2 and uJi > —d/2 for 
sufficiently large ^. These results are not consistent with 
the assumptions made in the derivation sketched above. 
Nevertheless, the zero- and slow-mode expansions seem 
to hold for all d > 2 and < ^ < 
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[40] Let us give this argument in more detail. The represen- 
tation B — SA in dimension d means that = djA^-' 
where A^^ — —A^^. The relation which replaces (|44p in 
general dimension d is 

where is the 2-point correlation of the 2-form A. 

The result that J d'^rC^-'{v) — follows if one assumes 
that the correlation function ^'^'-'^(r) — » sufficiently 
rapidly as |r| oo. The result (144 |l for d = 3 is recovered 
from the relation A^^ = e^'^^Ak between the 2-form A""^ 
and the usual vector potential Ak- 
[41] We remark that related eigenfunction expansions hold for 
the heat-kernels: 

j;^.^(r,0|p,t) = i?"(p,0|r,i) = ^e^°*f^*(r)fr,(p). 

a 

[42] Note that (57) applies for small r at fixed times t, whereas 
(56) applies at long times t for fixed r. However, the 
two results agree in their common domain of validity for 



r, L(0) < {DxtY^'' . This may be seen by applying (57) 
to W(i) to obtain for r < (-Dit)^/^ 

This result is verified in Section rill Cl with the explicit ex- 
pression for M^(i). Substituting the above into (|57p gives 

{r,t) ~ C(Dit)-(''+'^i+^2'/^Z('fj(r). 

This same result may be obtained by changing the inte- 
gration variable in (57) back to p —'p/ (Dit)^^'' and then 
employing the similar "quasi-equilibrium" result 

substituted into (57). 



